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Impact of CFO Estimation on the Performance 
of ZF Receiver in Massive MU-MIMO Systems 

Sudarshan Mukherjee, Saif Khan Mohammed and Indra Bhushan 

Abstract 

In this paper, we study the impact of carrier frequency offset (CFO) estimation/compensation on 
the information rate performance of the zero-forcing (ZF) receiver in the uplink of a multi-user massive 
multiple-input multiple-output (MIMO) system. Analysis of the derived closed-form expression of the 
per-user information rate reveals that with increasing number of BS antennas M, an 0{'/M) array gain 
is achievable, which is same as that achieved in the ideal zero CFO scenario. Also it is observed that 
compared to the ideal zero CFO case, the performance degradation in the presence of residual CFO 

(after CFO compensation) is the same for both ZF and MRC. 

Index Terms 

Massive MIMO, carrier frequency offset (CFO), multi-user, array gain, zero-forcing (ZF). 

I. Introduction 

In the past few years, massive multiple-input multiple-output (MIMO) systems have emerged 
as one of the key teehnologies in the evolution of the next generation 5G wireless systems due 
to their ability to support high data rate and improved energy efheieney [[H, lIH. In a massive 
multi-user (MU) MIMO system, the base station (BS) is provided with hundreds of antennas to 
simultaneously serve only a few tens of single-antenna user terminals (UTs) in the same time- 
frequeney resouree [l3l. Inereasing the number of BS antennas open up more available degrees 
of freedom, whieh helps aeeommodate more number of users, thus improving the aehievable 
speetral efheieney flUl, At the same time, the required radiated power to aehieve a hxed 
desired information rate ean be redueed with increasing number of BS antennas, M (array gain). 
It has been shown that even with imperfect channel state information (CSI), the achievable 
array gain for any sub-optimal linear receiver (e.g. zero-forcing (ZF), maximum ratio combining 
(MRC) etc.) is 0{VM) flU. 

Above results assume perfect frequency synchronization at the BS receiver, without which the 
performance of the system would deteriorate rapidly. In practice acquiring perfect knowledge 
of the carrier frequency offsets (CFOs) between the received user signals at the BS and the 
frequency of the BS oscillator is however a challenging task. There exists various techniques 
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for CFO estimation and compensation for conventional small MIMO systems in the literature 
na-Eoi. However these algorithms incur tremendous increase in computational complexity with 
increasing number of BS antennas, M and increasing number of UTs, K (i.e. massive MIMO 
scenario). Recently in IfTTI an approximation to the joint ML (Maximum Likelihood) CFO 
estimation has been proposed for massive MIMO system. However this technique requires a 
multi-dimensional grid search and therefore has high complexity with large number of UTs. 

In [fT^ the authors propose a simple low complexity algorithm for CFO estimation and a 
corresponding communication strategy for massive MU-MIMO uplink. It has been shown that 
with sufficiently large M, the algorithm has only 0{M) complexity (independent of the number 
of UTs). However the impact of the residual CFO (due to CFO compensation) on the performance 
of massive MIMO is yet to be studied. The most common linear suboptimal receivers used in 
massive MIMO uplink are MRC (maximum ratio combining) and ZF (zero-forcing) receivers. 
With the MRC receiver, system performance is limited by the multi-user interference (MUI) in 
the high SNR regime. For the ideal zero CFO scenario, the ZF receiver is known to remove this 
limitation by eliminating the MUI O . In this work we therefore study the impact of the residual 
CFO error (due to the CFO estimation strategy proposed in [fT^ l on the achievable information 
rate of the ZF receiver and compare it to that of the MRC receiver. To the best of our knowledge, 
this paper is the first to report such a study. 

The contributions of our paper are as follows: (i) we have derived a closed-form expression 
for an achievable information rate for the ZF receiver with MMSE (minimum mean square error) 
channel estimation and CFO compensation. A closed-form expression for the same is also derived 
for MRC; (ii) analysis of the ZF information rate expression reveals that an 0{y/M) array gain is 
achievable. This is very interesting since even for the ideal zero CFO scenario, the best possible 
array gain is known to be 0{'/M) only dH; (iii) for the same desired per-user information rate, 
the SNR gap (i.e. the extra SNR required by MRC when compared to ZF) does not degrade with 
CFO estimation/compensation, when compared to the ideal zero CFO case. This suggests that 
compared to the ideal zero CFO case, the performance degradation in the presence of residual 
CFO (due to compensation) is the same for both ZF and MRC. [Notations: C denotes the set 
of complex numbers. E denotes the expectation operator. (.)^ denotes the complex conjugate 
transpose operation, while (.)* denotes the complex conjugate operator. Also, 1^ denotes the 
N X N identity matrix and A^k (or (^)mfc) denotes the (m, fc)-th element of matrix A.] 

H. System Model 

We consider a frequency-flat massive MU-MIMO uplink (UL) channel, where the massive 
MIMO BS is equipped with M BS antennas and is coherently communicating with K single 
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CFO Compensation Strategy-I 

Data Communication Phase 


CFO Compensation Strategy-ll 

BS sends CFO estimates to the UTs for compensation 



Fig. 1 The communication strategy: CFO Estimation and Compensation Strategies and Data Communication. 

antenna UTs simultaneously in the same time-frequency resource. Therefore for a massive MU- 
MIMO BS, acquisition and compensation of CFOs from different UTs is important. Since in 
massive MIMO, the BS is expected to operate in time division duplexed (TDD) mode, the 
coherence interval (of Nc channel uses) consists of a UL slot (iV„ channel uses), followed by a 
downlink (DL) slot (Nc — Nu channel uses). As for the communication strategy (see in Fig. [I]), 
we perform CFO estimation in a special UL slot prior to the data communication. For CFO 
estimation, we adopt the CFO estimation strategy presented in ifT^ . CFO compensation can 
be performed in two different ways - (i) at the BS (prior to channel estimation and multi¬ 
user detection); or (ii) at the respective UTs prior to data transmission (this however requires 
transmission of CFO estimates from the BS to the UTs over a control channel in the DL slot, 
following the special UL slot for CFO estimation). Data communication starts from the first 
UL slot, following the special CFO estimation UL slot. In these UL slots, prior to UL data 
transmission, the UTs transmit pilots for channel estimation^ The special UL/DL slot for CFO 
estimation might be repeated every few coherence intervals, depending on how fast the CFOs 
change. 

A. CFO Estimation Strategy in ([72]/ 

For the CFO estimation phase, special pilots are transmitted by the UTs in the uplink. A pilot 
sequence of length N < Nu is divided into B = N/K pilot-blocks, where each pilot-block is 
K channel uses longl^ Each UT transmits only a single impulse of amplitude \/Kpa in each 
pilot-block. Therefore in the 6* pilot-block, the fc* UT transmits impulse at f = T{b, k) = 
(h — 1)K +k — l, where /c = 1, 2,..., iC and 6 = 1,2,...,!?. The pilot signal received at the m* 
BS antenna at time r(6, k) is therefore given by rm[r(6, k)] = Qmk +Wm[T(b, k)], 

where cok = 27rA/fcTs is the CFO for the /c* user (Tg = 1/Be, where Be is the communication 
bandwidth and A/^ is the frequency offset of the UT). Also gmk = hmk\/l\, m = 1, 2,..., M 
and k = 1,2,... ,K, is independent complex baseband frequency-flat channel gain coefficient 

'The same CFO estimates can be used for CFO compensation prior to precoding in the DL slot of each coherence interval. 

^Although the CFO estimation method assumes N/K to be integer, we can accommodate non-integer values of N/K, by 
defining the number of blocks B = (N/Ky Hence for non-integer N/K, the S"' block is less than K channel uses long. 
Therefore the effective N/K is \N/K'\ for UTs allowed to transmit in the B"' block, and it is {\N/K'\ — 1) for all other UTs. 
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between the m* BS antenna and the UT and hmk ~ CN{0, ^)u\fWk > 0 models the geometric 
attenuation factor for the user and Wm[T{h^ k)] ~ CAA(0, a^) is the complex circular symmetric 
AGWN noise with variance The estimate of the CFO of the k^ UT, is obtained as the 
principal argument of the block-wise correlation term of the pilot sequence received from the 

B-l M 

n n EE C[r(&,fc)]?'m['r(5-hl,fc)] 

fc* user, i.e., u)k = ^ arg {pk)n whereld pk = —ttt-;-= Gke^‘^>=^ + vk, and 


M 


MK{B - l)pu/ 


G. = -^7 E \hr.t\ 

M rn=l 


Remark 1. Note that the above CFO estimate is well-defined iff \uJkK\ < tt. For most practical 
massive MIMO systems, this condition will hold true ifT^ . Also, from the strong law of large 
numbers it can be shown that for i.i.d. hmk, ^ 1 as M —)■ oo. □ 

Result 1. {Approximation of the CFO Estimate in If \uJkK\ -C vr and 7 = > 7°, 

then the above CFO estimate can be approximated hy Uk = arg {pk) Uk + where 

7 ° =-p-yand = '^{uk). Note that {uk —ojk) ~ -^^(0, cr'^J, where 

AG,.[y/l + 2MgE^-lJ 

the mean square error (MSB) given by 

o-L =IE[(wfc 


al is 


M{N - K)K^Gl 


( 1 ) 


Remark 2. Clearly with M —)■ oo, we have 7° oc Therefore we choose some constant cq > 0 


such that 7 = :7> 7^^ as M —)■ 00, thereby satisfying the required condition 7 :§> 7*^ in Result 

1. From O we note that with 7 = and M —)■ 00 (fixed N, K), we have E[(wfe - rs 


lA 


2K^{N K)^" since ^im Gk = 1. This shows that the MSB for CFO estimation approaches a 

constant value as M —)■ 00 with 7 oc 1/\/M for fixed K and fixed N. Note that with M —)■ 00, 
the desired MSE oc I/cq, i.e., a smaller desired MSE can be attained using a higher value of cq. 
Therefore for a target/desired MSE, a sufficiently large M must be chosen so that the required 
power (oc 1/\/M) for CEO estimation is within the desired limits. □ 

B. Uplink Data Communication 


After CEO estimation, CEO compensation can be performed in one of the following two ways: 

(a) the BS can feed the individual CEO estimates back to the corresponding UTs over a control 

channel in the DE slot, following the special UE slot (see Eig. [U). In this way, the UT would 
correct its CEO by rotating the transmit signal at the channel use by However in this 

method there is possibility of corruption of the estimates due to error in the control channel; 

(b) another way of correcting frequency offsets is to perform CEO compensation at the BS, 

^Independent and identically distributed Rayleigh fading is a commonly used model for the distribution of channel gains in 
a massive MIMO system o-di. 

''Here arg (c) denotes the ‘principal argument’ of the complex number c. 

^For expression of Vk see IT2I . 
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prior to channel estimation and multi-user detection. In this paper, we would use this seeond 
technique for CFO eompensation and study the impaet of the residual CFOs (i.e. Qk — ^k) on 
the performance of massive MIMO uplink^ The uplink data communieation starts at t = 0 
(see Fig. [T]). We assume that in the first K eonseeutive ehannel uses, the UTs transmit pilots 
for ehannel estimation sequentially in timeO, i.e., the UT transmits an impulse of amplitude 
y/Kpu only in the {k — 1)* ehannel use. The reeeived pilot at the m* BS antenna at t = k — 1 
is therefore given by r^lk — 1] = y/Kpugmk + Wm[k — 1], where k = 1,2,..., K and 

m = 1, 2,..., M. Wm[k — 1] ~ CAA(0, cr^) is the eireular symmetrie AWGN. 

C. MMSE Channel Estimation 


To estimate the ehannel gains, firstly, CFO eompensation is performed on the reeeived pilots 
at the BS. Sinee the pilots from different UTs are separated in time, pilot from the fc* UT 
after eompensation is given by ym[k — 1] = ^mlk — = y/Kgmk + nmk[k — 1], 

where! gmk = gmk ~ CA/'(0,1), nmk[k -l]=Wm[k- ^ CN'{D, cr^) and 

Acjfc = cDfc — cjfc is the residual CFO error. Next we eompute the minimum mean square estimate 


(MMSE) of the effective ehannel gain coefficient pmk as given below: 

y/TCiEPk 


9mk — 


r2/m[fc - 1] = 


y/TCp^Pk 


+n,„fe[fc- 1]). 


( 2 ) 


KpuPk + cr'^ Kp,,Pk + 0-2 

where m = 1, 2,..., M and k = 1,2,..., K. Using (O, the estimate of the effective ehannel 


gain matrix is given by ^ ,- ~ 

G = {y/KiuG^^ + N)D, 

where G = \gmk\MxK, G = [gmk]MxK, and $o = diag{l,e-^^^^,- 
N = [nmk[k - 1 ]]mxk and D = {y/KpjK + where D = diag{Pi, ^ 2 , ■ ■ ■ ,Pk)- 

III. Uplink Receiver Processing 


(3) 

Here 


In this section we formulate a generalized approach towards multi-user receiver proeessing 
at the massive MIMO BS. From Fig. [T] it is elear that uplink data transmission begins at the 
t = channel use and eontinues till the — 1)'*’ ehannel use, where is the duration 
of the UL slot. Let y/p^Xk[t] be the information symbol transmitted by the A:* UT at the 
ehannel use. The signal reeeived at the m* BS antenna in the ehannel use is given by 

K 

Pm[f] = y/p^ gmq Xq[t\ + Wm[t], where 171 = 1,2,..., M. To detect information symbols 

9=1 

®Note that the information theoretic performance is identical for both the CFO compensation techniques. 

^Though impulse type pilots are not amenable to practical implementation (due to high peak-to-average-power ratio (PAR)), 
we use them because our main objective is to study the first order effects of system parameters M, K, N, p„ etc. on the 
information rate performance of ZF and MRC detectors in the presence of residual CFOs. 

*Both g^nk and Wm\k — 1] have uniform phase distribution (i.e. circular symmetric) and are independent of each other. 

Clearly, rotating these random variables by fixed angles (for a given realization of CFOs and its estimates) would not change 
the distribution of their phases and they will remain independent. Therefore the distribution of Ijmk and rimk [k — 1] would be 
same as that of gmk and Wm [k — 1] respectively. 
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of the UT, first, CFO compensation is performed, followed by detection using the detector for 
the UT. In this paper, we only consider linear detectors. We also assume ~ CAA(0,1), 
for f = iC, • • • , — 1 and are i.i.d. Let r[t] = (ri[f], r 2 [f], • • • ,rM[t])'^- The detected signal 

from the /c'*’ UT at the BS after CFO compensation is given by 


M 


K 


M 


Xk[t] = a)! r[t] e 


^-jUkt _ 


Detection .. 

compensation 


m—1 


—j^kt 


= V^E E <^*rnk9m, X, [i] 

< 7=1 \ m—1 / 


K 


q=i,q¥=k 


^ H 
= °’k9q 


f xS] + af nfc[t] 


M 

'y ] ^mk '9‘mk [t] 

m—1 

=af Tifc[t] 


K 


= Qke 

A 


H~ -jAujkit-(k-l)) 


Xk 


[t] + y/K E ^k9qe Xq[t] e 


j{aq-UJk)t 




(4) 


where = (ai^, 02 ^, • • • , aMkY ^ is the linear detector for the /c* user, Acj^ = ujk — (^k 

and nmk[t] = Also, nk[t] = {nik[t],n 2 k[t], ■ ■ ■ ,nMk[t]V, 9q = i9iq,92q, ■ ■ ■ ,9MqV 

(the g'*’ column of G) and Qg = (gig, g 2 q, ■■■ , gMqV = 9q 
A. Coding Strategy 

We define the effective channel estimation error as emk = dmk — gmk (see Section II-C). Let 
Cfc = (eifc, e 2 k, ■■■ , eMkV = 9k - gk, where gk = (gik,g 2 k, ■ ■ ■ ,gMkV column of G). The 
mean vector and the covariance matrix of are respectively given by E[efc] = 0 and 

y/KPu. 




Kpudk 
KpuPk + cr'^ 


- 1 


Kpudk + 0-2 


nk[k - 1] I (^gk - Qk'j 


H 


(a) 


dk cr 
Kpndk + cr"^ 


IM, 


(5) 


where (a) follows from Q and nk[k — 1] = {nik[k — 1], • • • , nMk[k — 1])^. Using gk = gk — ^k 
in we get 

Xk[t] = Xk[t] +a^nk[t] 


+ 


= E 


Vik^i {°'k{9q- ^q)^ 


= Sfc[t] 
K 


= ENfc[t] 




\q=l, 

q¥=k 


w) 


,j(u}q-€ik)t _ ^^^-jAujk(t-(k-l)) 


Xk[t] 


= MUIfc[t] 


Sk [t]] Xk [t] + (Sfe [t] - E [Sk [d] ) Xk [t] +MUIfc [d + ENfe [t ], 


( 6 ) 


= ESfc[t] 


= SIFfe[t] 


where SIFfc[l:] is the time-varying self-interfering component of the desired signal and SIFfc[f] + 
MUIfc[d -f ENfc[f] = Wk[t] is the overall effective noise term. Further E Sk[t] is the average 
value of S'fc[f], across several uplink data transmission blocks, i.e., several channel realizations, 
and is a function of t. The same is also true for the variance of 114[l]- Furthermore for a given 
t, across multiple uplink data transmission blocks, the realizations of Wk[t] are i.i.d. Hence 
for each channel use, t = K, K + 1,..., Nu — I, we have a additive noise SISO (single-input 
single-output) channel in ®. Thus for each user there are — K different SISO channels 
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with distinct channel statistics. Therefore we eonsider N^ — K ehannel eodes for eaeh user, one 
for eaeh SISO ehannel. The data reeeived for eaeh user in the ehannel use aeross multiple 
eoherenee intervals is jointly deeoded at the reeeiver [[T^ . This eoding strategy albeit not praetieal 
is useful in eomputing an aehievable information rateh 

IV. Achievable Ineormation Rate 

In essenee, from the above eoding strategy, we have — K parallel ehannel deeoders for 
eaeh user. For the f'*’ SISO ehannel of the user, we note that the eorrelation between the 
desired signal term ESk[t] and the overall effeetive noise Wk[t] is zero, i.e., from ® we have 


E 


ESl[t]Wk[t] 


E 




E 


\xk[t]\^[Sk[t] - E [Sk[t]\ } +xl[t]M\Jh[t] + xndENfc[t] 






E 


L 




= 0, since Xk [i] and Sk [i] are independent 

= 1,fc \ / 


= 0, since rik [i] is zero mean and' 
independent of Xk [i] 


xe- 


j{Qg-€)k)t 






— 0, since gfcand e^are orthogonal due to ' 
MMSE estimation and is function of 


= 0. (7) 


= 0, since Xi [f] are all i.i.d. 

where (a) and (b) follow from the definitions of SIFfc[f] and MUIfc[t] in ®. With Gaussian 
information symbols Xk[t], a lower bound on the information rate of the effeetive ehannel in 
® is obtained by eonsidering the worst case uncorrelated additive noise (in terms of mutual 
information), having the same varianee as With Gaussian information symbols, this worst 

ease uneorrelated noise is also Gaussian [fT4l . The varianee of Wk\i\ is given by E[|II4[f]p] = 
E[|SIFfc[f] +MUIfc[f] +ENfc[f]p]. Sinee all Xk[t] and nfc[f] are independent and zero mean, it ean 
be shown that E[SIF^[f]ENfc[f]] = E[MUI^[f]ENfc[f]] = 0. Also due to MMSE ehannel estimate it 
ean be shown that E[SIE^[f]MUIfc[f]] = 0. Therefore E[|lTfc 
E[|ENfc[f]p]. Also, E \ESk[t] 

rate is therefore given by the following lower bound on I{xk[t]; Xk[t]) 

I{xk[t];xk[t]) > log 2 (l + SINRfc[f]), where SINRfc[f] = E |ESfc[f]| 

th ..OO,- K, 10 


ESk[t] 

= E 

SlEk[t] 

= E 


= E 

1 1 
m 


2] =E[|SIEfc[f]|2]+E[|MUIfc[f]|2] + 

= 0. An aehievable 


E 




( 8 ) 


and the overall information rate for the user is thus given bwf 

= — log2(l + SINR4t]). 

t=K 


N„, 


(9) 


A. Mutual Information Analysis for the ZF Receiver 


Eor a ZE reeeiver, the deteetor matrix is defined as A = (ai, a 2 , • • • , a^) = G{G^G) 
Clearly, for ZE reeeiver, A^G = Ik, i.e., a^Qq = ^k,q = I if k = q and 0 if A; 7 ^ g, 

®In practice, coding could be performed across a group of consecutive channel uses within each transmission block, since the 
statistics of Wk [f] and Sk [f] would not change significantly within a small group of consecutive channel uses. 

*°In a wireless channel of bandwidth 200 KHz and a coherence interval of duration 1 millisecond, even with A = 10 UTs, 

the channel estimation overhead is only 5%. Further, CFO estimation is performed at a 5 to 10 times slower rate than channel 
estimation and therefore its overhead is expected to be less than 1% 1151 . We have therefore neglected the CFO estimation 
overhead in since it is a mere scaling factor, which does not impact the main conclusions of our work. 



































where q = 1,2,K and k = 1,2,... ,K. Substituting this result in ®, we get ESk[t] = 
Xk[t], where we have used the faet that 

Clearly, E[|ESfc[t]p] = Similarly E[|SIFfc[t]p] = , 

E[|MUI,[t]|2] =p„E \{{G^G)-^}kk\ E and E[|ENfc[t]|2] = a^E \{{G^G)-^}kk 


i=l 


Lemma 1. With MMSE ehannel estimates, it ean be shown that E 


{{G^G)-^} 


kk 


Pk 


+ 


'KpM 


/{M — K), where G is the MMSE estimate of effeetive ehannel gain matrix (see dH)). 


Proof: See Appendix A. ■ 

Proposition 1. Eor the ZE reeeiver, the lower bound in dS]) is given by the aehievable information 

rate -Rf [f] = log 2 (l + SINRf [f]), where 
SlNRf[i] = - 




1 — e 


-'ylAt-(k-iW 


A 


1 


K 

£if7ft + l ' 7, 


( 10 ) 


for E 

|MUIfc[f]|2 

and E 

|ENfc[f]p 

(see paragraph before 

E 

|ESfc[f]p 

, E 

\SlEk[t]? 

,E 

|MUIfc[f]p 

and E 

|ENfc[f]p 


M-K V/3fe ' 
where 7 = and f = A', A' + 1,..., — 1. 

Proof: Using the expression for E {{G^G)~^}kk from Lemma 1, we get the expressions 

Lemma 1). Using expressions of 
(see paragraph before Lemma 1) 
in the expression of SINRfc[f] in dH]) we obtain (fTO . ■ 

B. Mutual Information Analysis for the MRC Receiver 

For MRC receiver, A = G, or, Vfc = 1,2, ...,A". Substituting this 

result in dS), we get Sk[t] = ^/pi^\\gk\\'^ The desired signal ESfc[f] is 

therefore given by ESfe[d = E[5'fc[d]a;fcW = or, ESfc[d = 


Therefore E[|ESfc[t]|2] = p, (e [(G^G)^^]) 




Similarly, E[|SIEfe[t]|2] = 


E 


(G^G) 


kk 


(e [(G"G)fefe])^ 





K [- 

(G^G)fe, 

2-1 

Pu 

E E 



i—l,i^k ^ 




E 

2=1 


Kpu0i-\-cr‘^ 


E 


(G^G) 


kk 


, E[|MUIfe[t]|2] = 


and E[|ENfe[dP] = a2E (G^G)fcfc 


Lemma 2. With MMSE channel estimate G of the channel gain matrix G (see dH)), it can 

2l 


be shown that E 

2 


E 


G^G 


ki 


G G)kk 
= M 


= E 


Kp.0l 


Kp,if}k+o-'^ 


G^G 


kk 


= M{M + 1) 


KpuPk+a'^ 


and. 


KpuPk+o-^ I V irpuft+o-2 


Proof: See Appendix B. ■ 

Proposition 2. For the MRC receiver, the lower bound in dH]) is given by the achievable 

information rate R™'^[f\ = log 2 (l + SINR™'‘^[f]), where t = A', A' + 1,...,A„ — 1, and 


siNRrid = 


,-<At-{k-AY 


1 — e 




1 f 1 


1 


M \Pk KPil 


K 

E 


1 

7J 


( 11 ) 








































































and 


Proof: Firstly we substitute the expressions of E (G^G) 


kk 


E 


(G^G) 


ki 


E 


(G"G) 


kk 


and E 


from Lemma [21 in the expressions for E |ESfc [f] p 


,E 


|SIFfc[f]p 


,E 


IMUL 


|ENfc[f]p (see paragraph before Lemma 2 ). Using these in ([ 8 ]), we get (ITT]) . ■ 

Theorem 1. {Achievable Array Gain) Consider \uJkK\ <C tt, a fixed K, N (length of pilot 
sequence) and a fixed desired information rate for the f'*’ channel code of the fc* user 
and defined in Proposition [T| and Proposition [2] respectively). For both the ZF and MRC 

receivers, as M ^ oo, the minimum required SNR 7 to achieve the fixed desired information 
rate decreases as Alternatively, with M —)■ 00 and 7 oc 1/\/M, the achievable information 
rate for the channel code, i.e., R^f] or approaches a constant value. 

Proof: We have observed from Remark 2 that as M ^ 00 with 7 = -^= (constant cq > 0), 

A = 


the MSE for CEO estimation converges to a constant limiting value, i.e., lim 


= 


^/m 


Co > 0 (constant). Substituting this result in the expression for SINR^^[f] in (ITOl) and also in the 
expression for SINR™’'‘^ [t] in (ITT]) with 7 = we have 


lim SINR^'"[t] = lim SINR^*[d = 


=-Co(t-(fc-i))^ 


1 _ e-Co(i-(fe-i))^ 


> 0 (constant). 


( 12 ) 


Kbico 


2^2 


Erom ([T2l) it is clear that = log 2 (l + SINRf[f]) and Rf'^‘^[t] = log 2 (l + SINR™’'‘^[f]) 

would also approach constant limiting values as M —)■ cx) with 7 oc ■ 


Remark 3. Erom Theorem 1, it is clear that with every doubling in the number of BS antennas, the 
minimum required SNR to achieve a fixed per-user information rate decreases by approximately 
1.5 dB as long as the number of BS antennas M is sufficiently large. This shows that with the 
CEO estimation technique proposed in ifT^ . the ZE receiver (also the MRC receiver) yields an 
0{'/M) array gain in the massive MIMO uplink. This is interesting since even for the ideal 
zero CEO scenario with ZE/ MRC receiver, the maximum achievable array gain is known to be 
only 0{VM) d. □ 

Remark 4. Erom ([T2l) we have lim SINR™‘'[d = lim SINRffi]. Clearly, as M —)■ 00 with 7 oc 

M—^oo M—^co 

1 / '/M, the achievable information rate for both the ZE and MRC receiver approach the same 
lower bound. This shows us the new result that even with CEO estimation/compensation, MRC 
and ZE receivers have the same performance when M is sufficiently large. □ 

V. Numerical Results and Discussions 

We present a comparative discussion on the performance of the ZE and MRC receivers, 
with CEO estimation/compensation in frequency-flat massive MIMO uplink. Eor monte-carlo 
simulations, we assume an operating carrier frequency /c = 2 GHz and a maximum CEO of 1 
PPM of fc- The communication bandwidth is Be = 200 KHz. The coherence interval and the 
maximum delay spread are 1 ms and 5 /iS respectively. Thus \ijJk\ < ^ and Nc = \ms/Be = 200 
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channel uses. The duration of uplink is iV„ = 100 channel uses. The length of pilot sequence 
for CFO estimation is taken as N = 100 and the number of UTs is K = 10. At the start of 
each CFO estimation phase ujk assumes a random value uniformly distributed in Also 

for simplicity, we assume (3k = 1, Vfc = 1, 2,..., J-f. The information rate for each user is also 
computed analytically using Propositions 1 and 2 in (|9l) with = ]E[(cDfc — replaced by 
its approximation in ([T]) with = 1 (see Remark 1 ). 

In Fig. |2] we plot the variation of the minimum required SNR 7 = pu/cr^ (both analytical 
and simulated) to achieve a fixed information rate of 2 bpcu (bits per channel use) for the P' 
user versus the number of BS antennas, M (fixed iC = 10 and N = 100). Observe that the 
analytical approximation to the required 7 for both ZF and MRC is quite tight. Also for M 
sufficiently large, with every doubling in M, the required 7 decreases roughly by 1.5 dB (note 
the decrease in required SNR 7 from M = 320 to M = 640). This supports Theorem 1 and 
shows that with the discussed CFO estimation/compensation technique, an 0{y/M) array gain 
is achievable. Also note that the required 7 for ZF and MRC is the same for sufficiently large 
M > 320 (see Remark 4). However for finite M, ZF is more power efficient compared to MRC. 
For example at M = 80, ZF requires approximately 1.7 dB less power than MRC. Next we 
consider this extra SNR required by MRC when compared to ZF (denoted as SNR gap) for the 
same desired information rate Ji = 1, 2, 2.5 bpcu (for the P‘ UT) for fixed M = 80 and iC = 10 
(see Table HI)- From Table IH we make an interesting observation that the SNR gap between the 
ZF and MRC receivers is almost the same irrespective of whether we have the ideal zero CFO 


TABLE I SNR GAP BETWEEN ZF AND MRC Receiver for fixed M = 80, K = 10. 


Desired Per-User Information Rate 

SNR gap for Ideal/zero CFO case 

SNR gap with CFO compensation 

1 bpcu 

0.1 dB 

0.12 dB 

2 bpcu 

1.7 dB 

1.71 dB 

2.5 bpcu 

4.57 dB 

4.59 dB 
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scenario or the residual CFO (after CFO eompensation) scenario. Therefore the new result in 
this paper is that with CFO eompensation, there is no signifieant degradation in the SNR gap 
when eompared to the SNR gap in the ideal/zero CFO seenario. 

VI. Conclusion 


In this paper, we study the impact of low-complexity CFO estimation and compensation on the 

performanee of ZF receiver in massive MIMO uplink in a flat fading environment and compare 

it to that of the MRC reeeiver. The tight elosed-form analytieal expressions for information rates 

of ZF and MRC reveal that an 0{'/M) array gain is indeed achievable with CFO estimation. 

This is interesting since the best possible array gain for ideal zero CFO scenario is also known 

to be Finally the study of the SNR gap between ZF and MRC reeeivers for the same 

per-user information rate suggests that eompared to the ideal zero CFO ease, the performanee 

degradation due to residual CFO is same for both the ZF and MRC reeeivers. 

Appendix A 

Proof of Lemma 1 

From the relation Qmk = h^kV^ we have G = where H = [hmk]MxK- 

Clearly from © we have G = {^/K^G^^o + N)D = {y/K^H += 

= X =v 

+ V^XD = ZXD, where (a) follows from the faet that = Ik 


and D = {y/K^Ix + D Let rik and Vk be the fc* columns of N and V 

respectively. Sinee rNj = 1, 2, • • • , iC are all i.i.d. random veetors, 

Vk = D~^^'^)kk nk, fc = 1, 2, • • • , iC, are also independently distributed as 


iM) 


We also note that the eolumns of H and V are independent of eaeh other. Clearly, the 
same is also true for the eolumns of Z. Therefore we ean write Z = UQ, where Q = 
{KpuIk + and U = (-Ui, U 2 , ■ ■ ■ , uk), where Uk ~ CA/'(0, Im) Vfc = 1, 2,..., JT are 

i.i.d. random veetors. Now, we have | {G^G)~^\ = {{D^q D^/‘^QU^UQD^k^^QD)~^}kk, 

I J kk 

whieh follows from the faet that G = ZXD, X = and Z = UQ. Since D, Q 

\w-^)kk, where T = QD^/^^qD and 




-1 


and D are all diagonal, \ {G^G) * = \J.kk\ 

I J kk L 

w = u^u Wm{M, Im) is a K X K eentral Wishart matrix with M degrees of freedom. 

KpuPk + 


Clearly, E 

LI J kk- 




E[(W-i)fefe] = 


Apu/3fc + \ 


Apu/3: 


K 


triyv-^) 


(b) 


(M - K)Kp,Pl 


2 ’ 


where (6) follows from E[tr{W ^)] = EH. 


Appendix B 
Proof of Lemma 2 

From Appendix A we know that W = U^U is a eentral Wishart matrix with M degrees of 
freedom, i.e., Wkk is x^(2M) (ehi-squared) distributed. Therefore from definition of G we have 
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E 


E 


(G"G)fcfc ^^\nkrE[Wkk] = M 


KpnB 


KpuPk + 0-2 


E 


(G^g) 


ki 


^ E 




= m 


kk 


'E 


W, 


(G"G)fcfc 

2 “ 


ki 


(c) 


= M(M+ 1) 

KpnPl 


Kpuf3^ 


KpuPk + cr^ J \ Kp^Pi + ct 2 


KpuPk + 0-2 

KpuPi 


and 


M, 


where (a) and (6) follow from the faets that G = ZXD, X = Z = UQ and T = 

(see Appendix A). Also (c) follows from the faet that E[| W^p] = E[\uj^Ui\‘^] = M, 
sinee V/c = 1, 2 ,..., A' are i.i.d. CN'{0, Im)- 
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